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ABSTRACT 

We investigate the thermodynamic properties of D-brane— anti-D-brane pairs 
and non-BPS D-branes on the basis of boundary string field theory. We calcu- 
late the finite temperature effective potential of N D-brane— anti-D-brane pairs 
in a non-compact background and in a toroidal background. In the non-compact 
background phase transition occurs slightly below the Hagedorn temper- 

ature, and the D9-D9 pairs become stable. Moreover, the total energy at the 
critical temperature is a decreasing function of N as long as the 't Hooft cou- 
pling is very small. This leads to the conclusion that a large number N of D9-D9 
pairs are created simultaneously near the Hagedorn temperature. In the toroidal 
background case (M^g-D x Tp), a phase transition occurs only if the Dp-Dp pair 
is extended in all the non-compact directions, as long as the 't Hooft coupling 
is very small. The total energy at the critical temperature also decreases as N 
increases. We also calculate the finite temperature effective potential of non-BPS 
D-branes, and we obtain similar results. Then, we consider the thermodynamic 
balance between open strings on these branes and closed strings in the bulk in 
the ideal gas approximation, and conclude that the total energy is dominated by 
the open strings. 
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1 Introduction 



Understanding the properties of unstable D-brane systems, such as coincident D-brane— anti- 
D-brane pairs and non-BPS D-branes [1], has been a subject of much interest (for a review 
see, e.g., Ref. [2]). Type IIB string theory contains D-brane— anti-D-brane pairs of odd 
dimension and non-BPS D-branes of even dimension, whereas type IIA string theory contains 
D-brane— anti-D-brane pairs of even dimension and non-BPS D-branes of odd dimension. 
The spectrum of open strings on these unstable branes contains a tachyon field T. In such 
a brane configuration, we have T = 0, and the potential of this tachyon field has a local 
maximum at T = 0. If we assume that the tachyon potential has a non-trivial minimum, 
it is hypothesized that the tachyon falls into it. Sen conjectured that the tensions of these 
branes and the negative potential energy of the tachyon exactly cancel at the potential 
minimum [3]. This implies that these unstable brane systems disappear at the end of the 
tachyon condensation. 

If tachyon condensation occurs in a topologically non-trivial way, there remain some 
topological defects, such as kinks and vortices. We can identify these topological defects as 
lower-dimensional D-branes and the topological charge as the Ramond-Ramond charge of 
the resulting D-branes [1] [3]. These D-brane charges can be classified using K-theory [4] 
[5]. In particular, we can realize all the D-branes through tachyon condensation from the 
spacetime-filling branes, such as D9-D9 pairs and non-BPS D9-branes. 

In recent years, significant effort has been devoted to studying the time development 
of the tachyon field as a "rolling tachyon" [6]. The creation of closed strings in the rolling 
tachyon background has also been investigated recently [7] . Thus it is natural to consider the 
problem of how to realize the initial conditions of a rolling tachyon. In other words, we wish 
to consider how the tachyon field should be placed on the top of the potential. If there is a 
situation in which T = is first the stable potential minimum, namely, the D-brane— anti- 
D-brane pair or the non-BPS D-brane is stable, we can realize such a configuration. Then, if 
it becomes unstable, the tachyon will roll down to the potential minimum. We have pointed 
out in the previous papers that there are cases that these unstable branes become stable 
at very high temperatures [8] [9]. In particular, the spacetime-filling D9-D9 pair becomes 
stable at sufficiently high temperature in all cases we have studied. 

We have considered finite temperature D-brane— anti-D-brane systems [8] [9] 2 in the 
framework of boundary string field theory (BSFT) [10] [11]. This theory is based on the 
Batalin-Vilkoviski formalism [14]. For a superstring, the solution of the classical master 
equation is given by [15] [16] [17] 

Seff = Zq, (1.1) 

where S e ff is the spacetime effective action and Z Q is the disk partition function of the 
two-dimensional world sheet theory. If we calculate this solution in a constant tachyon 
background in the case of a single Dp-Dp pair in type II string theory, we obtain the tachyon 
potential [15] [17] [19] 

V(T) =2r p V p exp(-8|T| 2 ). (1.2) 
2 For related discussions see Refs. [12] and [13]. 
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Here, T is a complex scalar tachyon field, V p is the p- dimensional volume of the system that 
we are considering, and r p is the tension of a Dp-brane, which is defined by 

T P = (f) » ( L3 ) 

where g s is the coupling constant of strings. The minimum of this potential is at |T| = oo, and 
it satisfies Sen's conjecture [3]. In order to calculate the free energy by using the Matsubara 
method in the ideal gas approximation, we must calculate the one-loop amplitude. We have 
assumed that the relation (1.1) also holds if we include the higher loop correction of the 
world sheet theory, namely, 

Seff = Z, (1.4) 

where Z is partition function to all orders of the two-dimensional world sheet theory. If 
we consider the one-loop amplitude based on BSFT, we are confronted with the problem 
of the choice of the Weyl factors [20] [21] [22] [23], as we will explain in the next section. 
Andreev and Oft have proposed one choice [20], and we have computed the one- loop free 
energy using it [8] [9]. Then, we computed the finite temperature effective potential of D- 
brane— anti-D-brane systems near the Hagedorn temperature and discovered that there are 
cases in which T = becomes a stable minimum near the Hagedorn temperature, although 
it is a local maximum at zero temperature. The purpose of this paper is to generalize the 
above-mentioned calculation to finite temperature systems with multiple D-brane— anti-D- 
brane pairs and multiple non-BPS branes and to discuss their thermodynamic properties. 

This paper is organized as follows. In §2 we review the thermodynamic properties 
of strings near the Hagedorn temperature and those of open strings on the D-brane— anti- 
D-brane system [8] [9]. In §3 we evaluate the finite temperature effective potential of iV 
D-brane— anti-D-brane pairs in a non-compact background and in a toroidal background 
(Mi ]9 _ D x T D ) on the basis of BSFT. In §4 we generalize this D-brane— anti-D-brane case 
to the non-BPS D-brane case. Then, in §5 we study the thermodynamic balance of open 
strings on these branes and closed strings in the bulk. In §6 we consider the choice of the 
Weyl factors by investigating the case of another choice as an example. Finally, §7 contains 
conclusions and discussions. 



2 Thermodynamics of a String Gas near the Hagedorn 
Temperature 

In this section, we review the thermodynamic properties of an ideal gas of strings near the 
Hagedorn temperature. We summarize the basic properties of a string gas in §2.1, and we 
summarize the results of previous papers [8] [9] with regard to the thermodynamics of a 
D-brane— anti-D-brane system in §2.2. 
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2.1 Thermodynamics of a String Gas 



In the canonical ensemble method, all the statistical variables are derived from the partition 
function. The partition function Z(j3) can be obtained from the free energy F(j3) as 

Z(/3) = exp[-/3F(/3)]. (2.1) 

The free energy can be calculated using the Matsubara method; that is, the free energy is 
given by the path integral of connected graphs of strings on the space where the Euclidean 
time direction is compactified with a circumference equal to the inverse temperature, f3. 
Thus, it seems that we can calculate the statistical variables of strings by using the canon- 
ical ensemble method. However, we cannot trust the canonical ensemble method near the 
Hagedorn temperature for the following reason [24]. The partition function Z((3) is given by 
the Laplace transformation of the density of states fl(E): 

Z(P) = / dE Q{E)e- pE . (2.2) 
Jo 

We can expect that the canonical ensemble method gives the same values as the microcanon- 
ical ensemble method if the integrand has a sharp peak. However, the density of states of 
strings behaves as 

n(E) ~ e^ E (2.3) 
for large energy E, where (3 H is the inverse of the Hagedorn temperature, 

H = 2nV2^>, (2.4) 

and the integrand has no sharp peak near the Hagedorn temperature. For this reason, we 
must compute the finite temperature effective potential by using the microcanonical ensemble 
method, which is more fundamental than the canonical ensemble method in the sense that 
it is derived directly from ergodic theory. The energy E is more useful than f3 to analyze the 
string gas at high temperature, because f3 remains almost equal to f3 H , varying little even if 
E varies considerably. 

All the statistical variables are derived from the density of states in the microcanoni- 
cal ensemble method. The density of states, Q(E), can be obtained from the inverse Laplace 
transformation of Z((3) as 

/■L+ioo AR 

n(E) = Z z ^ e > ^ 

JL-ioc AlTl 

where we must treat (3 as a complex variable in this case. We must choose the constant 
L such that the path of integration in the complex /3-plane lies to the right side of all the 
singularities of the integrand. Then, we can deform the contour to the left so as to pick up 
the singularities, as sketched in Figure 1. The density of states has been derived from this 
formula in the case of closed strings [25] [26], in the case of open strings on D-branes [27], 
and in the case of open strings on a single Dp-Dp pair [8] [9]. 

In order to calculate the density of states from (2.5), we must investigate the singular 
part of the free energy, because the partition function is obtained from the free energy as 
(2.1). Let us separate the free energy into regular and singular parts, as 

F{f3) = F reg {f3) + F smg {f3) ) (2.6) 
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Figure 1: Complex (3 plane. 



and expand the regular part in (/3 — (3 H ) as 

-!3F reg ~ A V P - c7 V p (/3 - (3 H ). (2.7) 

It is easy to see that Ao is a constant with dimensions of inverse volume, and <jq is a constant 
with dimensions of energy density. We call <To the Hagedorn energy density. Then, we obtain 
the density of states as 

d(3 
lc a 2ni 



n(r, E) ~ e fe^+A v P /■ ^ exp r (/3 _ Ph)E + Fs . ng{(3) 



:2.8) 



where E = E — a V p and the contour C a is taken as sketched in Figure 1. A general property 
of the singular part of the free energy of strings is that its leading singularity is at (3 = (3 H . 
We call this the Hagedorn singularity. If we obtain the explicit form of Q(E) from (2.5), the 
entropy S(E) can be derived from Q(E) as 

S(E) =lntt(E)5E, (2.9) 

where 5E represents the energy fluctuation, and the inverse temperature (3 is given by the 
partial derivative of S(E) with respect to E: 

f=m _ (2 ' 10) 

The finite temperature effective potential of the Dp-Dp system can be calculated from these 
variables as 

V(T, E) = V(T) - i3~ 1 S. (2.11) 

From this, we can study the stability of a system at finite temperature. As we see below, 
there are cases in which unstable brane configuration at zero temperature becomes stable at 
finite temperature. 
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2.2 Finite Temperature Effective Potential for Single D-brane— anti- 
D-brane Pair 

In this subsection we review the results of our previous papers [8] [9]. As explained in the 
previous subsection, we must first compute the free energy of open strings in order to compute 
the finite temperature effective potential. It can be obtained by calculating the amplitude in 
the space where the Euclidean time direction is compactified with a circumference equal to 
the inverse temperature (3. Let us employ the weak coupling approximation and treat strings 
as an ideal gas. We take into account only the one-loop amplitude, with the corresponding 
cylinder world sheets winding in the Euclidean time direction at least once. 

The one-loop amplitude of open strings in a Dp-Dp system in the Minkowski space- 
time has been investigated on the basis of BSFT [20] [21] [22] [23]. The one-loop world sheet 
of an open string has the annulus topology and two boundaries. When we compute the 
one-loop amplitude of open strings in a Dp-Dp system, we are confronted with an ambiguity 
in the choice of the Weyl factors of the two boundaries of this world sheet. This is because 
the conformal invariance is broken by the boundary terms in BSFT. At the tree level with a 
disk world sheet, we do not have such a problem, because it has only one boundary. 

Andreev and Oft have proposed the following form of the one-loop amplitude [20] in 
the Minkowski spacetime in type II string theory on the basis of the principle that its low 
energy part should coincide with that of the tachyon field model [15] [17] [19]. If we restrict 
ourselves to a constant tachyon field, which we denote by T, it is given by 



where V p is the volume of the system that we are considering. This amplitude can be obtained 
straightforwardly by choosing the boundary action as 



Let us call this the 'cylinder boundary action'. This action is natural in the sense that both 
sides of the cylinder world sheet are treated on an equal footing in this case. Here we are 
considering only the coincident Dp-Dp system, and we do not treat the case of a parallel 
brane and antibrane pair with a finite distance [13] [28] in this paper. 

If we compute the one-loop amplitude in the space where the Euclidean time direction 
is compactified with a circumference equal to the inverse temperature /3, we obtain the free 
energy. The result is [8] 




(2.12) 




(2.13) 
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Figure 2: Here, for clarity, we display the brane and antibrane as if they were separated but, 
in fact, we consider the case of a coincident D-brane— anti-D-brane pair. 
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(2.14) 



This is the free energy in a non-compact flat background. This free energy can be obtained 
from the proper time form of the free energy [29], which is given by 



+ 



(2yra 
V, 



r 2 (3 2 



/ -(4vrr)-^ £ £(-1)^ -2™' 



iirar ) 



(2-KCt' 



by multiplying by 2 and by substituting the mass spectra 



M R 2 r- 



r 2 (3 2 



M, 



NS 



-[N B + N NS + 2T 
a' V 

±(n b + N r + 2\T\ 

rv' \ 



(2.15) 

(2.16) 
(2.17) 



where M^s and Mr are the mass of the Neveu-Schwarz and Ramond sectors, respectively, 
and Nb, Nns and Nr are the oscillation modes of the boson, Neveu-Schwarz fermion and 
Ramond fermion, respectively. This can be interpreted within open string theory as follows. 
As depicted in Figure 2, there are four types of open strings, namely, two types of strings 
whose two ends attach to the same brane and two types of strings that are stretched between 
the Dp-brane and the Dp-brane. If we denote the open strings extending from a Dp-brane to 
a Dp-brane as (Dp, Dp), the others are denoted as (Dp, Dp), (Dp, Dp) and (Dp, Dp). We must 
impose the GSO projection for (Dp, Dp) and (Dp, Dp) and the opposite GSO projection for 
(Dp, Dp) and (Dp, Dp). Therefore, instead of employing the GSO projection, we may multiply 
by an overall factor of 2. 

The mass spectra given in (2.16) and (2.17) are shifted by 2\T\ 2 /a f from the mass 
spectra on the Dp-Dp system at T = 0. This shift of mass is reminiscent of the mini- 
superspace approximation in S-brane thermodynamics proposed by Maloney, Strominger 
and Yin [30]. Their argument is based on boundary conformal field theory (BCFT), and the 
boundary action is chosen as that corresponding to the case of a time dependent tachyon 
field. 3 They treat it as a time dependent mass and compute the creation rate of the open 
strings. Our model is related to their model in the sense that the boundary term is treated 
as a shift of the mass spectra. 

We can also compute the free energy in a flat background with some toroidally 
compactified directions in a similar way. We assume that D-dimensional space is compactified 
and that the rest of the (9 — _D)-dimensional space is left uncompactified (namely, M lj9 _ D x 
Tp). We call this background the 'toroidal background'. We also assume that the Dp-Dp 
system extends in the d- dimensional space in the non-compact direction and in the (p — d)- 
dimensional space in the toroidal direction. The Dp-Dp system extends in all non-compact 
direction if D + d — 9. The free energy is given by 



F(T,P,R) = 



16tt 4 V, 




A 



2iR/T 



a' 



(2.18) 



where Vd is the ci-dimensional volume in the non-compact directions parallel to the Dp-Dp 
system. The effect of compactification is reflected by the infinite product of ^-functions. 
We can also obtain this free energy from (2.15) by multiplying by 2 and substituting the 

3 The boundary action in the bosonic string case is chosen as 

Sb = A / dr cosh — ^J- , 



■I 



where A is a constant real number [30]. 
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mass spectra 



M, 2 



M, 2 



p—d 

R 2 



where mj denotes the momentum number in the directions parallel to the Dp-Dp system, 
and rii denotes the winding number in the directions transverse to it. It should be noted 
that these free energy and mass spectra are invariant under the T-duality transformation 
[31] 



a' 



Ri^^ , rnj^m, {2.21) 



for directions parallel to the Dp-Dp system and 



a' 



Ri^-E- , rn^mj, (2.22) 
til 

for directions transverse to it. From these free energies, we can compute the finite tempera- 
ture effective potential by using the method presented in the previous subsection. 



We now summarize the results of previous papers [8] [9]. In the non-compact back- 
ground case, the results for the D9-D9 case differ from those in the other cases. In the D9-D9 
case, the sign of the coefficient of the |T| 2 term of the finite temperature effective potential 
changes from negative to positive slightly below the Hagedorn temperature as the temper- 
ature increases. This means that a phase transition occurs near the Hagedorn temperature 
and the D9-D9 pair becomes stable. On the other hand, the coefficient remains negative in 
the Dp-Dp case with p < 8, and thus a phase transition does not occur. This leads us to 
the conclusion that only a D9-D9 pair is created near the Hagedorn temperature. In the 
toroidal background case, a phase transition occurs near the Hagedorn temperature only in 
the D + d = 9 case, namely, the case in which the Dp-Dp pair is extended in all non-compact 
directions. The coefficient remains negative in the Dp-Dp case with D + d < 9, and thus a 
phase transition does not occur. 



3 Finite Temperature Effective Potential for Multiple 
D-brane— anti-D-brane Pairs 



We have considered only a single D-brane— anti-D-brane pair to this point. Now we generalize 
the calculation of the previous section to the case of multiple D-brane— anti-D-brane pairs. 
We begin by reviewing the tachyon field in the case of multiple D-brane— anti-D-brane pairs 
in §3.1. We then investigate the behavior of the free energy near the Hagedorn singularity 
in §3.2. Finally, we calculate the finite temperature effective potential in a non-compact 
background and in a toroidal background in §3.3 and §3.4, respectively. 
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3.1 Tachyon Field for Multiple D-brane— anti-D-brane Pairs 



In the single D-brane— anti-D-brane case, the tachyonic field is a complex scalar field. If we 
consider open strings on N D-brane— anti-D-brane pairs, the spectrum of the open strings 
contains a tachyon field in the (N, N) representation of the U(N) x U(N) gauge group [32]. 
In this case, the tree level tachyon potential in BSFT can be written [15] [17] 



V(T) = 2t p V p Trexp (-8T f T) , 



(3.1) 



where T is an N x N complex matrix. This potential depends on T only in the form T^T. 
At the one-loop level, the free energy can be written 

F(T, P) = T^T Trexp (-2vrTtTr) Tr exp (-2^^) 

Ph j T 2 



X 



' M0\ir) 



i(3 2 



if3 2 



(3.2) 



This free energy also depends on T only in the form T^T. Because T^T is a hermitian 
matrix, it can be diagonalized by a unitary transformation. This corresponds to choosing 
the appropriate linear combination of branes and antibranes. Let us choose T^T as 



/ in 



v 



o 



\ 



|T| 2 J 



(3.3) 



which corresponds to the choice that all the tachyon fields on each Dp-Dp pair condenses by 
the same amount. The finite temperature effective potential does not depend on the choice 
of T, as long as it gives above form of T^T. Thus we can choose T as 



/ T 



\ 



(3.4) 



V T J 

Substituting (3.4) into (3.1) and (3.2), we obtain 

V(T) = 2iVr p V p exp(-8|T| 2 ) 

and 



(3.5) 
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l6ir 4 N 2 V., 



F{T,0) 



A 



p+i 



-/ — exp(-47r|T| 2 r 

JO T ~ V 



X 




- 1 



(3.6) 



respectively. The tachyon potential (3.5) is N times that of a single Dp-Dp pair, and the 
free energy (3.6) is N 2 times that of a single Dp-Dp pair. This free energy can be obtained 
from the proper time form of the free energy (2.15) by multiplying by 2N 2 and substituting 
the mass spectra (2.16) and (2.17). The factor of 2N 2 comes from the fact that there are 
2N 2 types of strings with (Dp, Dp) and (Dp, Dp) and 2N 2 types of strings with (Dp, Dp) and 
(Dp, Dp). The 't Hooft coupling g s N must be very small for the ideal gas approximation in 
this case. 4 



It should be noted that we can reproduce the tachyon potential and free energy on 
(N — n) pairs with arbitrary positive integer n by changing n of the diagonal components of 
the matrix (3.4) to oo. Thus, we may consider a very large matrix from the beginning, and 
we can treat all iV cases in the same framework by making such choice of the matrix. 



3.2 Free Energy near the Hagedorn Singularity 



As mentioned in §2.1, the leading singularity of the free energy is at P = f3 H) which is 
called the Hagedorn singularity. This can be seen from the free energy (3.6) by changing the 
variable of integration from r to t, given by 



1 



(3.7) 



and considering the region of large t. Then, using the modular transformation of i? functions, 
we obtain 



F(T,0) 



Pu 



p+l 



r°° . p-9 / 
J at t 2 exp I 

'# 3 (0|it)V 



47r|Tr 



t 



( MOW ] 
W(0\it)) 



#3 



if3 2 t 



Pi. 



#4 I 



ip 2 t' 



P, 



- 1 



(3.8) 



4 If the 't Hooft coupling is very large, we must consider the black brane solution in supergravity, which 
corresponds to the brane-antibrane system [18]. However, because we are assuming that the 't Hooft coupling 
is very small, the perturbative calculation in open string theory should be valid. 
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Expanding the d functions and extracting the leading term in the large t region near the 
Hagedorn singularity, we obtain 



4N 2 V 



PS 



p-9 



wr A dtt 2 ex P 



4tt|T| 

T 



IT- 



P 2 ~P„ 



2 1 



Pn 



(3.9) 



where we have introduced the low energy cutoff A. This cutoff is required, because we cannot 
apply the above approximation to the small t region. Because we will compute the |T| 2 term 
of the finite temperature effective potential in the vicinity of T = 0, we expand the free 
energy in |T| 2 and keep the lower order terms: 



F(T,/3) 



AN 2 V P r°° ,. .p^9 



Pn 



■ I l . , , ■ ( P Ph , 
^pf / dtt 2 exp -7r — — t 



P» 



16nN 2 V p \T\ 2 r°o ^ / p 2 -p 



P 



p+i 



dt t 2 exp — 7r- 



P, 



(3.10) 



We can rewrite (3.10) near (3 = P„ in terms of the incomplete V function as 



P. 

2 a+5 7r a+2 N 2 V p (P - p H ) a+1 \T\ 



A 



9-a 



P, 

r ( -a- 1 , 2vr 



A 



(3.11) 



where we have defined 



_ 7-p 
a = — - — . 



(3.12) 



We can deduce the singular part of the free energy F sing from (3.11) as follows. 



(a) p = 9 (a = -1) 

When p — 9, which means that a — — 1, the first argument of the incomplete T function 
in the first term of (3.11) becomes I. Therefore we can set A = 0. For the second term 
of (3.11), we can use the formula 

T(0, x) = - 7 - Ins - E Hrd"' (3 - 13) 
where 7 is the Euler constant. Combining these two terms, we get 

w IT m 2N ' V - 16nN 2 V r \T\ 2 , / 13-Jh \ 

F -'" (T - i}) = — ln r ~ A J ■ (3 - 14) 

To make our notation the same as that of Ref. [8] we set A = (27r) _1 in this formula 
hereafter. 



12 



(b) p : even (a : half- integer) 



When p is even, the first arguments of the incomplete T functions are negative half- 
integers. In this case, we can set A = and we obtain the singular part of the free 
energy as 

Ph 



8(27r) a+2 r(-a - l)N 2 V p \T\ 



2 



(0 - (3.15) 



(c) j9 : odd (a : integer) 

When p is odd, the first arguments of the incomplete T functions are negative integers. 
The incomplete T function whose first argument is a negative integer can be expanded 

as 

1 n (—l) n 
r(-n, x) = —e~ x ^(-l) 5 - 1 ^ - s)\ x- n+s - 1 + ^-i-r(0, x). (3.16) 
n. s=1 n. 

In our case, A — > corresponds to x — > 0, and therefore the incomplete V function 
can be approximated by the last term. Thus, the singular part of the free energy is 
obtained as 



(3.17) 

We will also set A = (27r) _1 in this formula hereafter. 



From these, we can calculate the finite temperature effective potential by using the 
method described in §2.1. For N pairs of Dp-Dp, (2.7) is modified as 



-PF reg ~ A iV 2 V p - a N 2 V p (P - p H ), 



(3.18) 



and (2.8) as 



Q(T, E) ~ gfes+^v, ^ ^ exp _ + Fs . ng(/3) ] j (3 10) 



where we have defined the quantity E as 

E = E- a N 2 V p . 



(3.20) 



From these, we compute the finite temperature effective potential near the Hagedorn tem- 
perature. 
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3.3 Non-compact Background 



In this subsection, we compute the finite temperature effective potential of multiple D- 
brane— anti-D-brane pairs in a non-compact background. Although the method of calculation 
depends on p, we explicitly consider only three cases as examples, since calculation is almost 
the same as that in Ref. [8]. 



(a) p = 9 (a = -1) 



Let us first consider the case of N pairs of D9-D9, which is the most interesting case. 
In this case (3.19) can be rewritten as 



/3 h E+\ N 2 Vq 



Q(T,E) ~ e 



where we have defined 



dP (l3-f3 h 



-N 2 D-!V 9 \T\ 2 



c a 2iri 



exp 



P~P, 



(3.21) 



C-i 



_/3 8' 
167T 



P, 



(3.22) 
(3.23) 



Let us suppose that both E and E/V 9 are very large. Then, the saddle point method 
works well, because the exponent in the integrand is very large. The result is 



tt(T,E) 



1 / P H 2 E Y N2D - lV » m2 [N^VX* 



E° 



2^ \N 2 C_ L V 9 



x exp p H E + X N Vq + 2NJC_ 1 V 9 E 



(3.24) 



From this Q(T,E) and the equations (2.9) and (2.10), we can obtain the entropy 
S(T, E) and the inverse temperature P as 



S(T, E) 



P 



NtD.MlTl 2 ^ ( p H 2 E 



3 -m 



E 



{N^C^VgJ 4 \NtC^V 9 i(5E) 
+P H E + X N 2 V 9 + 2Ny/c_ 1 V 9 E, 



2l 



(3.25) 
(3.26) 



respectively. It is noteworthy that the temperature is lower than the Hagedorn tem- 
perature if the energy density is very large and T = 0. The finite temperature effective 
potential can be derived from (2.11). In order to demonstrate the stability of the D- 
brane— anti-D-brane system, we need only the |T| 2 term of V(T, E). This term is given 

87riY 2 V 9l (nP H w E\\ 2 



-16iVT 9 V9 + 



Ph 



10 



In 



2N 2 V 9 



(3.27) 
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It should be noted that the second term in the coefficient of |T| 2 is an increasing 
function of E. Because the first term is constant as long as Vg and r 9 are fixed, the 
sign of the |T| 2 term changes from negative to positive at large E. The coefficient 
vanishes when 

27V 2 V 9 f2p H w r 9 \ 



E ~ 



7T& 



TIN 



If we approximate (3.26) as 



(3.29) 



we can derive the critical temperature T c at which the coefficient vanishes. The result 
is 

-l 

1 + exp 



T c ~ AT 1 
= AT 1 



1 + exp 




(3.30) 



where we have used the explicit formulae (1.3) and (2.4) in the second equality. Here, 
we see that this temperature is very close to the Hagedorn temperature, since r 9 is 
very large if the coupling of strings g s is very small. It should be noted that the 
critical temperature decreases as iV increases. Above this temperature, the coefficient 
of |T| 2 is positive and T = becomes the potential minimum. This implies that a 
phase transition occurs at the temperature T c , which is slightly below the Hagedorn 
temperature, and the D9-D9 system is stable above this temperature. We can obtain 
the results for the single D9-D9 pair case considered in Ref. [8] by substituting N — 1. 

From (3.20) and (3.28), we can calculate the total energy at the critical temperature 

as 



E tot ~ a N 2 V 9 + 2Nr 9 V 9 + 



2N 2 V 9 



10 



exp 



'2(3 



H r 9 



TlN 



= a N 2 V 9 + 



2NVa 



+ 



2N 2 V 9 



(2tt)V 5 # s 7r(87r 2 a') 5 



■ exp 



(3.31) 



where we have included the brane tension energy, and we have also used the explicit 
formulae (1.3) and (2.4) in the second equality. It should be noted that for small 
values of A" and g s , the total energy at the critical temperature is a decreasing function 
of N. This implies that the multiple D9-D9 pairs are created simultaneously not in 
succession. Let us define N m i n as the value of A" that minimizes the total energy 
at the critical temperature. Then the N min pairs of D9-D9 are first created as the 
energy of the system increases. If we ignore the first term of (3.31) and compute N min 
numerically, we obtain 



N ■ ~ 

1 v mm — 



45.29 

9s 



(3.32) 



However, for the ideal gas approximation, we must impose the condition that the 't 
Hooft coupling is very small, namely, 



g s N < 1, 



(3.33) 
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and iV cannot be taken arbitrarily large, even if g s is very small. Equation (3.32) does 
not satisfy the condition (3.33). Therefore, we cannot determine the value of N min 
with a perturbative calculation. We can only say that a large number N of D9-D9 
pairs are created simultaneously, because the total energy at the critical temperature 
is a decreasing function of N as long as the 't Hooft coupling is very small. 

A similar argument has been made by Danielsson, Guijosa and Kruczenski using a 
calculation of the finite temperature effective potential in the case that there is only a 
tachyon field [12]. In that case, the finite temperature effective potential is a decreasing 
function of TV if P is fixed in the canonical ensemble method, and it seems that N 
increases without bound. However, they asserted that, given a finite amount of energy, 
we can determine the value of N that maximizes the entropy and minimizes the finite 
temperature effective potential if we employ the microcanonical ensemble method, 
because the entropy vanishes if we use the entire energy to create D-brane— anti-D- 
brane pairs or if we create no pairs at all. Our calculation demonstrates this with a 
calculation including the string massive mode. However, we know of no way to obtain 
the value of N under the condition (3.33). 



(b) p = 8 (a = -1/2) 



Next, let us consider the case of N D8-D8 pairs. Here, we find a result strikingly 
different from that in the case of the D9-D9 system. From (3.19), the density of states 
is given by 



n(T,E) 



D f3 H E+\ N 2 V s 



o (3„E+\ N 2 V 6 



d£ 

C a 27T2 



cxp 



(p - P H )E + N 2 C_ h V 8 (P - P H )^ 
-N 2 D_MT\\f3-f3 H y-~ 



dp 
c a 2iri 
x exp 



1-N 2 D _MT\ 2 (P-P H )* 



(P - P H )E + N 2 C i V 8 (P - P H )~ , (3.34) 



where we have used the small \T\ 2 approximation in the second equality, and we have 
defined 



C_ 



D_ 



2§ 



P, 



2if7r 2 
' Pji 



(3.35) 
(3.36) 



We can also use the saddle point method, which yields 



Sl(T,E)~ - ;\_ 8 5 exp\p H E + \ N 2 V 8 + " ^ ^ 



\ 



35237T2E 6 



2§ 



2\E' A 



(3.37) 

We can calculate the entropy S, the inverse temperature P and the potential V(T, E) 
from (2.9), (2.10) and (2.11), as in the D9-D9 case. The |T| 2 term of V(T, E) is given 



/ 



16 



by 

T| 2 . (3.38) 

It should be noted that the second term in the coefficient of \T\ 2 decreases as E becomes 
large. Thus, the coefficient of |T| 2 remains negative for large E. This implies that a 
phase transition does not occur, unlike in the D9-D9 case. For the p = 7, 6 (a = 0, 1/2) 
case, we can calculate the finite temperature effective potential in a similar way, again 
using the saddle point method. 



-16iVr 8 V 8 -J- 

38 H 12 E* 



p = 3,l {a = 2,3) 



In this case, we cannot use the saddle point method, because we cannot ignore the 
corrections coming from the 0(V p (/3 — f3 H ) 2 ) term in the expansion of the regular part 
of the free energy F reg (3.18) if we use the saddle point method, as explained in Ref. 
[8]. We therefore adapt another type of approximation as follows. 

The density of states can be obtained from (3.19) as 



Q(T,E) 



D (3 H E+\ Q N 2 V P 



dp 
c a 2m 



exp 



W-Ph)e 



+ (N 2 C a ([3 - (3 H ) a - N 2 D a \T\ 2 {(3 - (3 H ) a+1 ) V p ln 



where we have defined 



^+i4(27r) a 



C a — 



D n = 



r(a + i)A/- Q ' 

(-l) a+2 8(27r)»+ 2 
T(a + 2)P H 8 ' a ' 



Then, using the transformation 



z = -(0-p H )E, 



we get 



q(t,e) ~ 



,J3 H E+X N 2 V P 
1 



dz 



c b 2m 



: cxp 



— z 



+ [ {-iy N 2 C a ^ - {-l)^N 2 D a \T\ 2 ^ 



z 
E 



a+l 



The contour C b is taken as sketched in Figure 3, where 

Zl = -{L' -/3 H )E. 



(3.39) 

(3.40) 
(3.41) 

(3.42) 



(3.43) 



(3.44) 
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We can approximate this integral with the following contour integral on both sides of 
the cut: 



Z1 dz _ z 
— e exp 

7T 



y a+l 



x sin 



-TT [ (-l)«jV 2 C a ==^ - (-l) a+1 N 2 D a \T\ 2 ^ 



E 



a+l AT 2 , 



y a+l 



+1 



(3.45) 



Expanding the integrand here in powers oi 1/E and taking the limit zi — > oo, this 
integral becomes 

4(27r) a N 2 V P 8(27r) a+2 |T| 2 iV 2 V p 



ir(3 H 7 - a E 
Thus, the density of states is given by 

\a 7\r2 



— 4- 



7r(3 H 8 - a E a+1 ' 

8n 2 \Tr 



^ _M2nTN% ePHE+XoN , Vp f 
[3 H 7 - a E a+1 V 

Then, substituting il(T,E) into (2.9), we obtain 



(3.46) 



(3.47) 



8n 2 \T\ 2 



S(T, E) ~ -(a + 1) In ^- , + + A iV 2 V p - — r^f, (3.48) 



and from (2.10), we get 



P 



a + l 8n 2 \T\ 



(3.49) 



From this, we see that the temperature is higher than the Hagedorn temperature if 
the energy density is very large and T = 0. This property is similar to that of a closed 
string gas in a non-compact background, as we see in §5. From the above, we can 
calculate the finite temperature effective potential, and its \T\ 2 term is given by 



16tt 2 

■16Nt p V p + jj= 



I T I ^ 



(3.50) 



From this, we see that the coefficient remains negative for large E, and therefore no 
phase transition occurs. For the p = 5,4,2,0 (a = 1,3/2,5/2,7/2) cases, we can 
calculate the finite temperature effective potential in a similar way. 



For other cases, we only present the results, because these can be obtained with 
calculations similar to that given in Ref. [8]. The |T| 2 term of the finite temperature 
effective potential is given by 



/V £ \AN 2 V 7/ 



\T\ 



(3.51) 



for the D7-D7 case, 



-16iVr 6 V 6 



2 u n b N*V 6 
p H 28 E 3 



(3.52) 
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Figure 3: Complex z plane. 



for the D6-D6 case, 



for the D5-D5 case, where 



and 



-16iVr 5 V 5 



16tt 2 



\T\ 



'N 2 V, 



-16Nt p V p + 



16tt 2 



f3 H 2 E 



T 



(3.53) 



(3.54) 



(3.55) 



for D4-D4, D2-D2 and DO-DO cases. From these results, we can see that the coefficients 
remain negative for large E, and therefore no phase transition occurs. We can obtain the 
results for the single Dp-Dp pair case considered in Ref. [8] by substituting N — 1. 



3.4 Toroidal Background 



In this subsection, we generalize the calculation given in the previous subsection to the case 
of a toroidal background. We consider a Dp-Dp system in type II string theory compactified 
on a .D-dimensional torus Td and assume that the Dp-Dp system is extended in the d- 
dimensional non-compact directions and in the (p — d) -dimensional compact directions. In 
this case, the free energy is given by 
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F(T,/3,R) 



16n 4 N 2 V d dr 



p-d D 



Ph 



d+1 



f4^"n n M° 

J0 T 2 I=li=p-d+l \ 



X 



tf.s 



p-d+l 



i^ 2 



#3 



2iiL 2 r 



a' 



PhT. 



- 1 



04 



- 1 



(3.56) 



This free energy is iV 2 times that of the single pair of Dp-Dp given in (2.18). We can 
obtain the free energy near the Hagedorn singularity with a calculation similar to that in 
the previous subsection. The result is 



F(T,P,R) 



p-d D / 

xn n M° 

i=i i =p -d+i \ 



2a' 



Ph' 

ia't 



2Ri 2 



An\T\ 2 N 2 AV v , p+d-ii 



Ja 



Ph 

p-d D I 

*n n M° 



i=i i= p -d+i 



dt r 



%Rj 2 t 



where we have defined 



and 



a 



2a' 

id-p+% 



#3 



cxp y n - 

ia't 



P 2 ~ P* 



P, 



2Ri 



2- 2 p H a ur= P -d + iRi 

p—d 

v P = v d n Ri- 

i=i 



(3.57) 



(3.58) 



(3.59) 



This free energy is invariant under the T-duality transformation given in (2.21) and (2.22). 
We only need to investigate the region in which R > \fa' . 



Let us first consider the case that all the radii are close to the string scale, y a' . In 
this case, the Hagedorn singularity is dominant, and we need to consider only the terms 



F(T,P,R) 



N 2 AV n f 00 , D +d 



Ph 



dt t 2 exp — 7r 



P 2 ~ Ph' 



Ph' 



t 



+ 



4n\T\ 2 N 2 AV n r°° 



P* 



Ja 



D+d-U 



P I dtt — 

A 



exp — 7T- 



p 2 -p L 



Pb 



(3.60) 



This free energy is proportional to that in the non-compact background case, (3.10), if we 
replace (D + d) by p. From this, we can conjecture that a phase transition occurs only when 
D + d = 9. In the D + d = 9 case, the |T| 2 term of the finite temperature effective potential 
is given by 



-IQNt p V p + 



27iN 2 AVp ( 2nP H E Y 
P H n \N 2 AV p ) 



(3.61) 
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It should be noted that the second term in the coefficient of |T| 2 is an increasing function of 
E. Because the first term is constant as long as V p and r p are fixed, the sign of the \T\ 2 term 
changes from negative to positive at large E. This means that a phase transition occurs near 
the Hagedorn temperature. The coefficient vanishes when 



- N 2 AV P (8/3 H r p \ 
E ~ — : — tt- 1 - exp 1 1 



\itNA ) 



(3.62) 



The critical temperature T c and the total energy E tot at the critical temperature are given 
by 

' irNA j 



1 + exp 



(3.63) 



and 



E tot ~ a N 2 V p + 2Nr p V p + exp 



^ uNA ) ' 



(3.64) 



respectively. The total energy at the critical temperature is also a decreasing function of 
N as long as the 't Hooft coupling is very small, as in the case of D9-D9 in a non-compact 
background. This implies that a large number N of Dp-Dp pairs are created simultaneously. 



In the D + d < 8 case, the |T| 2 term in the finite temperature effective potential is 



given by 



for D + d = 



for D + d = 7, 



for D + d = 6, 



for D + d = 5, where 



and 



■16Nt p V p 



2 5 tt 2 nIaIv, 



-16Nt p V p 



\T\ 



8n 2 N 4 A 2 V n \ ( f3 H E 



\N 2 AV P 



\T\ 



-1QNt p V p + 



2 w ir 6 N s A 4 V p 4 
6f3 H 4 E 3 



16tt 2 

— lOT v V v H n 

p 1 f3 H E 1 



\T\ 



E^E- 2 ^ A ^ln(N 2 AV p ), 



-1QNt p V p + 



16tt 2 

J7b 



(3.65) 



(3.66) 



(3.67) 



(3.68) 



(3.69) 



(3.70) 



for D + d < 4. From these results we can see that the coefficients remain negative for large 
E, and therefore no phase transition occurs in these cases. We can obtain the results for the 
single Dp-Dp pair case considered in Ref. [9] by substituting N — 1. 



Next, let us consider the case that the radii Bq in the d- dimensional directions 
transverse to the Dp-Dp system are much larger than the string scale, \fa'. In this case, the 
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Figure 4: Radius dependent singularities. 



radius dependent singularities exist near the Hagedorn singularity, as sketched in Figure 4, 
and we must consider the radius dependent terms. Because one of the infinite product of 
radius dependent d functions in (3.57) can be rewritten as 



p~d+d / 

n m° 

i=p-d+l \ 



ia't\ ~ / ^fn,V, 
—y = ^ exp - ^ —j-t \ , (3.71) 

Zrl d / n i= -oo \ i=p-d+l Zri d 



the free energy can be approximated as 



N 2 AV r°° ( 
F(T,fi,R) *> dt J2 exp -7 

— ^ exp -7T t , (3.72) 

iA * n l= -oo \ Ph J 



4n\T\ 2 N 2 AV p r°°dt 




where we have defined 

/ p-d+d 2 A 

^n 2 ^^ 2 1- E ^ ■ (3-73) 

\ i=p-d+l ZJX d J 

From this we see that singularities exist at (3 — (3 n on the real axis of the complex /3-plane, 
and they approach the Hagedorn singularity (3 = (3 H as increases. 

Even if R^ is much larger than the string scale, a phase transition occurs only in the 
D + d = 9 case. The coefficient of the |T| 2 term of the finite temperature effective potential 
becomes positive only when the density of states is dominated by the contribution from the 
Hagedorn singularity and we can use the saddle point method to compute the density of 
states, as in the N = 1 case, which we considered in Ref. [9]. These conditions are realized 
if E satisfies the relation 



E ^> max 



I a' IL ^r A R- d N 2 R- 4 ~ d V 

13.71) 



Then the finite temperature effective potential has the same form as that in the case that 
d = and D + d = 9. A phase transition occurs when E is given by (3.62), which satisfies 
(3.74) as long as the 't Hooft coupling is very small. The critical temperature is given by 
(3.63). Because A depends on Rq, the total energy at the critical temperature is an increasing 
function of Rj as long as the 't Hooft coupling is very small. 
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4 Non-BPS D-brane 



We briefly study the non-BPS D-brane case in this section. If we consider N non-BPS D- 
branes, the spectrum of open strings contains a tachyon field in the adjoint representation 
of the U(N) gauge group [1]. The tension of a single non-BPS Dp-brane is given by V2t p 
[33], and the tree level tachyon potential of iV non-BPS D-branes is given by [15] [17] 



V(T) = v / 2iVr p V p exp(-2T 2 ), 



(4.1) 



where T is a real variable, and we have chosen the matrix as (3.4). This is l/y/2 of the 
tachyon potential in the Dp-Dp, case with |T| 2 replaced by T 2 /4. The one-loop free energy 
is given by half of that in the Dp-Dp, case with |T| 2 replaced by T 2 /4 [20]. The free energy 
of open strings on iV non-BPS D-branes, for example, in a non-compact background is given 
by 



F(T,P) = 



^+3-exp(-7rT 




(4.2) 



This free energy can be obtained from the proper time form of the free energy (2.15) by 
multiplying by N 2 and substituting the mass spectra (2.16) and (2.17) with 2|T| 2 replaced 
by T 2 /2. The factor N 2 can be interpreted in the context of open string theory as follows. 
There are N 2 types of strings projected by the ordinary GSO projection and iV 2 types of 
strings projected by the opposite GSO projection [1]. Therefore, instead of employing the 
GSO projection, we may multiply N 2 as an overall factor. From these considerations, we 
see that we can obtain the results in the case of non-BPS D-branes from those in the case of 
D-brane— anti-D-brane pairs by replacing N by Nj\f2. Below, we give the results only for 
the cases of non-BPS D9-branes in a non-compact background and of non-BPS Dp-branes 
with D + d = 9 in a toroidal background, in which a phase transition occurs. 

In the case of N non-BPS D9-branes in a non-compact background, the T 2 term of 
the finite temperature effective potential is given by 



« /om v , 47riV 2 Vg _ ( ir(3 H 10 E 
-8v^vr 9 V 9 + -^ln(^^ 



T 



(4.3) 



The critical temperature T c and the total energy E tot at the critical temperature are given 
by 5 



1 + exp 



ttN 



(4.4) 



5 We have not divided the first term on the right hand of (4.5) by y/2, because we have defined cr by 
(3.18) and E by (3.20). 
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and 



E tot ~ a N 2 V 9 + V2Nt 9 V 9 + 



10, 



N 2 V 9 / gyg^ -T,, 



■ cxp 



(4.5) 



_/3 10 

In the case of AT non-BPS Dp-branes with D + d = 9 in the toroidal background, the T 2 
term of the finite temperature effective potential is given by 



(4.6) 



The critical temperature T c and the total energy E tot at the critical temperature are given 
by 



1 + exp 



ttNA 



and 



E tot * . iV 2 V p + V2Nr p V p + ^ exp 



(4.7) 



(4.J 



In these two cases, the total energy at the critical temperature is also a decreasing function 
of N as long as the 't Hooft coupling is very small, as in the Dp-Dp case. This implies that 
a large number N of non-BPS D-branes are created simultaneously by a phase transition in 
these cases. No phase transition occurs in the other cases. 



5 Thermodynamic Balance 



To this point, we have ignored closed strings. We need to consider not only open strings 
but also closed strings if we wish to deal with the entire the system of strings in type II 
string theory. Abel, Barbon, Kogan and Rabinovici have investigated the finite temperature 
system of BPS D-branes on a 9-dimensional torus and treated the thermodynamic balance 
between open strings on the BPS D-branes and closed strings in the bulk [27]. In a similar 
way, in this section we treat the thermodynamic balance between open strings and closed 
strings in the case that D-brane— anti-D-brane pairs or non-BPS D-branes become stable 
near the Hagedorn temperature. 

We begin by reviewing the thermodynamics of closed strings in a 10-dimensional non- 
compact spacetime. We can reach the Hagedorn temperature for closed strings by supplying 
a finite amount of energy if we consider only closed strings [34]. Based on this fact, it has 
been said that the Hagedorn temperature in the closed string case is associated with a phase 
transition, in analogy to the deconfining transition in QCD. Sathiapalan [35], Kogan [36] and 
Atick and Witten [37] have argued that this phase transition occurs because the 'winding 
modes' of the Euclidean time direction become tachyonic above the Hagedorn temperature, 
and these tachyon fields condense towards the potential minimum. This phase transition is 
called the 'Hagedorn transition'. We do not yet know where the minimum of this tachyon 
potential is and in what kind of backgrounds the system becomes stable. Thus, there is a 
possibility that we cannot apply the ideal string gas approximation to closed strings above the 
Hagedorn energy density. However, we temporarily assume that we can treat closed strings 
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as an ideal gas, because it is expected that closed strings cannot reach the Hagedorn energy 
density if we consider thermodynamic balance between closed strings and open strings, as 
we now explain. 



We can investigate the thermodynamic balance of open strings and closed strings by 
considering the fact that the entropy becomes maximal in the equilibrium state. Let us fix 
the total energy E tot and consider the variation of the total entropy S tot . E tot is given by 

E tot = E + E c + E b , (5.1) 

where we have denoted the energy of open strings by E a , that of closed strings by E c , and 
that of the brane tension by E b . S tot is given by 

Stot — S + S c , (5-2) 

where we have denoted the entropy of open strings by S and that of closed strings by S c . 
By using the definition of the temperature (2.10) in the microcanonical ensemble framework, 
we can calculate the partial derivative of S to t with respect to E c as 

dStot = dS^_dS^_ _L 1_ 

f)P flP f)P T T ' v 0,0 / 

UHi c UHi c UHi 1 dosed J-open 

where we have used the fact that E b is a constant and the approximation that E tot is fixed in 
the first equality. When Stot realizes its maximum and the system becomes the equilibrium 
state, this partial derivative vanishes. Therefore, the thermodynamic balance condition is 
given by 

Topen Tclosed- (^ - ^) 

That is, thermodynamic balance is realized when the temperatures of the two types of strings 
become the same. If the temperature of one component is higher than that of other, then 
energy flows from the hotter to the cooler. 

Let us consider, for example, the D9-D9 system in a non-compact background under 
the condition that the energy is sufficiently large to create D9-D9 pairs. The entropy of open 
strings is given by (3.25), and it can be approximated for large E as 



S ~ (3 H E + 2N^/C_ 1 V 9 E (5.5) 
at T = 0, where E a = E Q — a N 2 V p . The temperature is given by 



T ~ 

■'■open — 




(5.6) 



which is lower than the Hagedorn temperature, i.e., Topen < T H . We need an infinite energy 
to reach the Hagedorn temperature, and this is a 'limiting temperature' [34] in the case of 
an open string on the D9-D9 pair. On the other hand, if we calculate the entropy of closed 
strings that have an energy density that is much larger than the Hagedorn energy density in 
the microcanonical ensemble method, we obtain [25] 6 

S c ^P H E c -^\n( fl] E \ \, (5.7) 

z \ Vg " 5E C ii 



6 Because the free energy of a closed string gas behaves as if it were the case that a — 5, we can calculate 
the entropy and the temperature as in case (c) of §3.3. 



25 



where E c ' = E c — ooVg and 5E C represents the energy fluctuation. The temperature is given 




which is higher than the Hagedorn temperature, i.e., %i OS ed > Th, and the temperature 
approaches the Hagedorn temperature from above. Therefore, the thermodynamic balance 
condition (5.4) is not satisfied as long as the closed strings have an energy density larger than 
the Hagedorn energy density, and the energy flows from the closed strings to the open strings. 
To put it another way, by substituting these temperatures into (5.3), we can conclude that 
S t ot is a monotonically decreasing function of E c throughout the entire Hagedorn regime. 
These facts imply that energy flows from closed strings to open strings before the closed 
strings realize the Hagedorn energy density. The closed strings cannot realize the Hagedorn 
energy density, and its temperature does not reach the Hagedorn temperature. Therefore, 
if we consider thermodynamic balance, open strings dominate the total energy of strings as 
long as this total energy is large enough to create D9-D9 pairs. 

If we consider the thermodynamic balance of open strings on the Dp-Dp pairs with 
D + d = 9 and closed strings in a toroidal background, we can also show that open strings 
dominate the total energy of strings [27], although there are cases in which the temperature 
of the closed strings is lower than the Hagedorn temperature even if they have an energy 
density larger than the Hagedorn energy density. This is because open strings in these cases 
have a larger entropy than closed strings if they have the same energy, and energy flows 
from closed strings to open strings in order to maximize the total entropy [27]. 7 We can 
derive a similar conclusion in the non-BPS D-brane case. Therefore, for all cases in which 
a phase transition occurs, open strings dominate the total energy, and closed strings have 
only a small energy. These facts indicate that the open string degrees of freedom are very 
important near the Hagedorn temperature. 



6 Annulus Boundary Action 

In this section, we consider an annulus-type world sheet as an example of a different choice 
of the Weyl factors. In the single Dp-Dp pair case, the natural choice of the Weyl factors of 
this world sheet gives the following form of the boundary action [21] : 8 



7 D9-branes and non-BPS D9-branes in a non-compact background and Dp-Dp pairs and non-BPS Dp- 
branes with D + d = 9 in a toroidal background are classified into L[— 1] in Ref. [27]. L[— 1] systems dominate 
the string energy, because they have the largest entropy under the condition that each type of string have 
the same energy. They are the most 'favorable'. 

8 We can choose another type of Weyl factor for the annulus world sheet [22] In this case, we obtain the 
results similar to those for the choice (6.1). 




(6.1) 
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where we have set the radius the inner boundary to q and that of the outer boundary to l. 9 
We call this action the 'annulus boundary action'. The one-loop partition function is given 
by [20] [21] 
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dq 
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If we change the variable of integration from q to r, given by 

7T S 



q = exp — 



(6.2) 



(6.3) 



we can obtain a form of the one-loop amplitude that is similar to that in the case of the 
cylinder boundary action (2.12). Explicitly we have 



_ 1GkHV p 

Zj\ — 



Pa 



p+1 



1 dr 
^rrexp 
o r — 



[1 + e 



'M0\itY 



' M0\ir) ' 
*i'(0|«-). 



(6.4) 



If we compute the one-loop free energy using the Matsubara method, we obtain 
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(6.5) 



We cannot obtain this free energy by substituting any form of the mass spectra into the 
proper time form of the free energy (2.15). Therefore we cannot explain this free energy in 
terms of a simple mass shift, like (2.16) and (2.17). In order to see the singular behavior of 
the free energy near the Hagedorn singularity, let us change the variable of integration from 
t to t, according to (3.7), and extract the leading term in large t region near the Hagedorn 
singularity. We then obtain 



F(T,/3) ~ 



4V, 



" e -m 2 



p+i 



A 



J at t 2 exp I — 7r — — — -i 



-a , 2n — A 

Pa i 



(6.6) 



where A is the low energy cutoff and a is defined by (3.12). It should be noted that the 
contribution of the tachyon T comes from only one boundary. This is very unnatural and 
a very different situation from that for the cylinder boundary action, which gives the same 
contribution to both boundaries. If we expand this free energy in |T| 2 , the power of (f3 — f3 H ) 

9 If we choose the radii as qr and r, with a positive real number r, respectively, we can absorb r through 
the redefinition of T. 
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in the \T\ 2 term is different from that in (3.11). Therefore, we obtain a different result 
than in the cylinder boundary action case. The change from the power in the cylinder 
boundary action case to that in the annulus boundary action case effectively corresponds to 
the replacement of a by (a — 1), so that a phase transition occurs in the Dp-Dp case with 
p < 7. We only present the results for N Dp-Dp pairs in these cases. In the D9-D9 case, the 
|T| 2 term of the finite temperature effective potential is given by 



-1Wt 9 V 9 + 



! N 2 VaE 



\T\ 



and the total energy at the critical temperature is 

E tot ~ a N 2 V 9 + 2Nt 9 V 9 + 2 9 vr^ 10 r 9 2 V 9 . 
In the D8-D8 case, the |T| 2 term is given by 

2*N*V a *E* 



-16Nt 8 V 8 + 



3& 



T 



and the total energy at the critical temperature is 
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In the D7-D7 case, the \T\ 2 term is given by 



4iV 2 V 7 / 3 H 8 E 
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and the total energy at the critical temperature is 
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(6.10) 



(6.11) 



(6.12) 



The total energy at the critical temperature is a increasing function of N in the D9-D9 case, 
while it is a decreasing function of iV in the D8-D8 and D7-D7 cases, as long as the 't Hooft 
coupling is very small, as in the cylinder boundary action case. However, we cannot obtain 
the value of N min for D8-D8 and D7-D7 pairs in the weak coupling region (3.33). Thus we 
cannot determine which Dp-Dp pairs are created first as the total energy increases. In the 
case of lower-dimensional Dp- Dp pairs, no phase transition occurs. In a toroidal background, 
a phase transition occurs only in the D + d > 7 case. We obtain a similar results in the 
non-BPS D-brane case. 



The property common to the cylinder boundary action case and the annulus bound- 
ary action case is that only the higher- dimensional branes are created by the phase transition 
in a non-compact background. Whether a phase transition occurs or not depends on the 
power of (f3 — f3 H ) of the |T| 2 term in the free energy near the Hagedorn singularity. It seems 
that this power in a non-compact background is a decreasing function of p for any choice of 
the Weyl factors. In the cylinder boundary action case, for example, this power is (a — 1), 
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as we can see from (3.11), and a is defined as (3.12), so that the power is a decreasing 
function of p. This property originates from the contribution from the momentum modes 
of open strings, which can be taken in the p-dimensional spatial directions on the Dp-Dp 
pairs. Therefore, it is expected that only higher- dimensional branes are created by the phase 
transition for any choice of the Weyl factors. It is also expected that only branes with large 
(D + d) are created by the phase transition in a toroidal background for any choice of the 
Weyl factors. 



7 Conclusion and Discussion 



In this paper, we have discussed the thermodynamic properties of Dp-Dp pairs and non-BPS 
D-branes in a constant tachyon background. We generalized the argument given in Refs. [8] 
and [9] to that for N Dp-Dp pairs and evaluated the |T| 2 term of the finite temperature 
effective potential by using the microcanonical ensemble method near the Hagedorn tem- 
perature. In the D9-D9 case in a non-compact background and in the D + d = 9 case in 
a toroidal background, a phase transition occurs slightly below the Hagedorn temperature. 
The total energy at the critical temperature is a decreasing function of N as long as the 
't Hooft coupling is very small in these cases. This implies that a large number N of D- 
brane— anti-D-brane pairs are created simultaneously. However, we cannot determine the 
value of N m i n with a perturbative calculation. In order to determine it, we must perform a 
non-perturbative calculation based on, for example, the matrix model [38] or the IIB matrix 
model [39]. The K-matrix model may also be useful, as it explicitly contains the tachyon field 
[40]. Such a phase transition does not occur in the other cases. We have also investigated 
the behavior of a finite temperature effective potential on non-BPS D-branes and obtained 
similar results. 

The following is another reason that we need to perform a non-perturbative calcu- 
lation. We have treated the phase transition near the Hagedorn temperature by calculating 
the |T| 2 term of the finite temperature effective potential up to the one- loop effect of open 
strings. By analogy with the Hagedorn transition in closed string theory, it seems that we 
cannot ignore the higher-loop contribution to the finite temperature effective potential if the 
total energy is larger than the brane tension energy. However, it is not clear whether we 
can ignore the higher-loop effect or not, because we are considering the |T| 2 term of the 
finite temperature effective potential of open strings. Actually, we can ignore the one-loop 
contribution in comparison with the tree level contribution in the lower- dimensional brane 
cases, as we have seen, for example, in §3.3. Thus, we must compute the higher-loop effect 
explicitly in order to verify our calculations up to one-loop. However, we have not succeeded 
in computing the higher-loop contribution by using BSFT to this time, because we do not 
know the natural choice of the Weyl factors in the higher-loop case, unlike in the one-loop 
case. We leave the non-perturbative calculation for the future. 

We have investigated the thermodynamic balance between the open strings on these 
branes and the closed strings in the bulk in the ideal gas approximation, and we have found 
that the open strings dominate the total energy. Thus, we conclude that the open string 
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degrees of freedom are very important near the Hagedorn temperature. Strictly speaking, 
however, we must consider the Boltzmann equation for open strings and closed strings in 
order to elucidate the thermodynamic balance, even if we can treat the system perturbatively 
[41]. 

The calculations given here are based on the cylinder boundary action, which was 
proposed by Andreev and Oft [20]. As mentioned in §2.1, there is a problem in choosing 
the Weyl factors in the two boundaries of the one-loop world sheet. The cylinder boundary 
action is natural because with it, both sides of the world sheet are treated on an equal footing, 
and its low energy part coincides with that of the tachyon field model [15] [17] [19]. If we 
calculate the finite temperature effective potential in the annulus boundary action case, we 
obtain results that differ from those in the cylinder boundary action case. However, in a non- 
compact background, a phase transition also occurs only in the case of higher-dimensional 
branes. As mentioned in §6, it is expected that only higher-dimensional branes in a non- 
compact background and branes with large (D + d) in a toroidal background are created by 
the phase transition, with any choice of the Weyl factors. 

It is possible that the D-brane— anti-D-brane pairs and the non-BPS Dp-branes are 
stable in the early universe, as the temperature is extremely high in this period. Then the 
universe expands in an inflationary manner, because the tension energy of these branes can 
provide an effective cosmo logical constant. A considerable number of studies have been made 
on brane inflation [42] [43] [44] (for a review see, e.g., Ref. [45]). We must consider gravity 
coupled with the tachyon field, including the finite temperature tachyon potential, in order 
to analyze a time dependent background. If we use the finite temperature effective potential 
we obtained, it is natural to treat this background in the framework of BSFT. The time 
dependent background at zero temperature was calculated by Sugimoto and Terashima in 
the framework of BSFT [46]. It would be interesting to generalize their calculation to the 
finite temperature case. It is expected that, even if these branes are stable initially, they 
become unstable, and the tachyon starts to roll down from T = [6], because the energy 
density decreases as the universe expands. 

It is noteworthy that spacetime-filling branes, such as the D9-D9 pairs and non-BPS 
D9-branes, are created at sufficiently high energy not only in a non-compact background but 
also in a toroidal background, since these branes always satisfy D + d = 9. These spacetime- 
filling branes are very advantageous in the sense that all the lower-dimensional D-branes 
are realized as topological defects through tachyon condensation from the spacetime-filling 
branes [4] [5]. It would be interesting to consider the possibility that our 'brane world' forms 
as a topological defect in the process of tachyon condensation in a cosmological context [44]. 

Finally, the phase transition to spacetime-filling branes is reminiscent of the phase 
transition in the Plank solid model of Schwarzschild black holes [47] . It might be interesting 
to study black holes as multiple spacetime-filling branes. 
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